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■ Coupled second order nonlinear differential equations are of fundamental impor- 
I tance in dynamics. In this part of our study on the integrability and linearization of 

■ nonlinear ordinary differential equations wc focus our attention on the method of 
deriving general solution of two coupled second order nonlinear ordinary differential 
equations through the extended Prelle-Singer procedure. We describe a procedure 
to obtain integrating factors and required number of integrals of motion so that the 
general solution follows straightforwardly from these integrals. Our method tackles 
both isotropic and non-isotropic cases in a systematic way. In addition to the above, 
we introduce a new method of transforming coupled second order nonlinear ODEs 
into uncoupled ones. We illustrate the theory with potentially important examples. 
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, 1. Introduction 

o 



In this part of our study on the integrability and linearization of nonlinear ordinary 
differential equations (ODEs) we focus our attention on the theoretical formulation 
and applications of the modified Prelle-Singer (PS) procedure (Prelle & Singer 1983; 
] Duarte et al. 2001; Chandrasekar et al. 2005a; 2006) to a set of two coupled second 

order ODEs. The need for this demonstration is due to the fact that classifying and 
studying two degrees of freedom dynamical systems are highly nontrivial problems 
' in the theory of nonlinear dynamical systems. Historically, several techniques have 

been proposed to identify and obtain general solutions of two coupled second order 
ODEs. To cite a few we mention Painleve analysis. Lie symmetry analysis, general- 
ized Noether symmetries technique, direct method and so on (Ramani et al. 1989; 
Lakshmanan & Sahadevan 1993; Bluman & Anco 2002; Lakshmanan & Rajasekar 
2003). Each of these methods have their own advantages and limitations. For exam- 
ple, among the above, certain methods fulfil necessary conditions alone whereas the 
others guarantee only sufficient conditions for the complete integrability of the sys- 
tem concerned. This factor alone is a motivating factor to search for more and more 
powerful methods to isolate and classify integrable and non-integrable dynamical 
systems. In this direction, by extending the PS procedure and its applications to 
coupled second order ODEs we argue that the PS method can be used as a stand- 
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alone technique to solve a wide class of ODEs of any order irrespective of whether 
it is a single or coupled equation. 

We mention here that the present analysis is not a straightforward extension 
of the scalar case. In fact by prolonging the theoretical formulation to the coupled 
second order ODEs, we deduce the determining equations for the integrating factors 
and null forms appropriately such that one can obtain the aforementioned functions 
in a more efficient and straightforward way. Thus the method of obtaining the 
integrating factors for the given equation is also augmented in this procedure in 
an efficient manner. Further, while studying the coupled dynamical systems one 
may face both isotropic and non-isotropic cases. Our method covers both of them 
in a natural way. In addition to the above, in this paper, we also introduce a new 
method to transform two coupled second order ODEs to two uncoupled second 
order ODEs. Thus, the PS procedure inherits several remarkable features both at 
the theoretical foundations as well as in the range of applications which we have 
listed out already in Chandrasekar et al. (2005a). Finally, we note that we have 
carefully fixed the examples so that the basic features associated with this method 
and the results which it leads to could be explained in an efficient way. 

The plan of the paper is as follows. In §2 we describe the PS method applicable 
for coupled second-order ODEs and indicate the new features in finding the inte- 
grating factors and integrals of motion. In §3 we establish a connection between 
integrating factors and the form of equations. In §4, the imcoupled equations are 
briefly considered. In §5, we discuss elaborately the method of constructing integrals 
and general solutions for the coupled nonlinear ODEs. We support the theory with 
two nontrivial examples which are being discussed in the contemporary literature. 
We also briefly discuss the application of our procedure for the case of Liouville 
integrable systems in §5d. We devote §6 to demonstrate yet another method to 
identify transformation variables from the first integrals which can be effectively 
used to rewrite the system of coupled ODEs into uncoupled ones so that one can 
integrate the resultant equation easily and obtain the general solution. We present 
our conclusions in §7. 



2. Prelle-Singer method for coupled second order ODEs 

(a) General Theory 
Let us consider a system of two coupled second order ODEs of the form 

^=^ = 1^' = i = l,2, (2.1) 

where Pi and Qi are analytic functions of the variables t,x,y,x and y. Let us 
suppose that the system (2.1) admits a first integral of the form I{t,x,y,x,y) = C 
with C constant on the solutions so that the total differential gives 

dl = Itdt + Ixdx + lydy + I±dx + lydy = 0, (2.2) 

where subscript denotes partial differentiation with respect to that variable. Rewrit- 
ing (2.1) in the form 

P P 

-^dt -dx = 0, -p^dt -dy = (2.3) 
Qi Q2 
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and adding null terms si(t, x, y, i, y)idt— Si{t, x, y, x, y)dx and S2{t, x,y,x, y)ydt — 
s^it, X, y, X, y)dy with the first equation in (2.3), and u\{t, x, y, x, y)xdt—ui{t, x, y, x, 
y) dx and U2{t,x,y,x,y)ydt — U2{t,x,y,x,y)dy with the second equation in (2.3), 
respectively, we obtain that, on the solutions, the l-forms 
p 

(— — h six + S2y)dt — sidx — S2dy — dx = 0, (2.4a) 
Qi 

P2 

{— — h uix + U2y)dt — uidx — U2dy — dy = 0. (2.4b) 
Q2 

At this stage, we wish to point out that one can also analyse the coupled second 
order ODEs (2.1) by rewriting them as a set of four coupled first order ODEs of 
the form x = xi, Xi — Pi/Qi, y = yi, yi = P2/Q2 and its equivalent one forms. 
By introducing four integrating factors, one can deduce the relevant determining 
equations by following the procedure given by us in the earlier paper, Part-III 
(Chandrasekar et al. 2008), to the above system of first order ODEs. However, 
after examining several examples we find that it is more advantageous to solve the 
system (2.1) in the second order form itself rather then introducing more number 
of variables. The procedure is as follows. 

Now, on the solutions, the l-forms (2.2) and (2.4) must be proportional. Mul- 
tiplying (2.4a) by the factor R{t,x,y,x,y) and (2.4b) by the factor K{t,x,y,x,y), 
which act as the integrating factors for (2.4a) and (2.4b), respectively, we have on 
the solutions that 

dl = R{(t)i + Sx)dt + K{<i)2 + Uy)dt - RSdx - KUdy - Rdx - Kdy = 0, (2.5) 

where 0i = Pi/Q^, i ^ 1,2, and S = ^"'+^"^ and U = ^"^+^"^ . Comparing 
equations (2.5) and (2.2) we have, on the solutions, the relations 

It = R{(j>i + Sx) + K{cl)2 + Uy), h = -RS, ly = -KU, h = -R, h = -K. (2.6) 

The compatibility conditions between the equations (2.6), namely Itx = Ixt, hx = 

^xt ? ^ty ^yt ) ^ty ^yt ? ^xy -^yx ) ^xx ^xx 7 ^yy ^yy 7 -^xy -^yx ) -^yx ^xy 

and Ixy = lyx, provide us the conditions, 

D[S] = -<l)ix - ^hx + ^S<P2x + ScP,^ + (2.7) 

D[U] = -K-^<f>iy + ^Uct>,y + U<f>2y + U^ (2.8) 

D[R] = -{R(f>ii + Kcf>2x + RS), (2.9) 

D[K] = -{K<j)2y + Rd^iy + KU), (2.10) 

SRy = -RSy + UK^ + KUx, Rx=SRi + RSi, (2.11) 

Ry = UKi + KUi, Kx=SRy + RSy, (2.12) 

Ky = UKy + KUy, Ry = K ^ . (2.13) 

Here the total differential operator, D, is defined by-D = ^-|-a;^-|- y-^ + M 
(^2^- Integrating equations (2.6), we obtain the integral of motion, 



/ = ri + r2 + rs + r4 



K + — { n + r2 + r3 + ri 



dy, (2.14) 
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where 



n =j(^R{cj>i+Sx)+K{cj>2 + Uy)jdt, r2 = - J (^RS + £{n)^dx 
ra =-j(KU+i-{n+r2)]dy, = - J 



dx. 



Solving the determining equations, (2.7)-(2.13), consistently we can obtain expres- 
sions for the functions {S,U,R,K). Substituting them into (2.14) and evaluating 
the integrals we can constract the associated integrals of motion. It is also clear 
that equation (2.1) can be considered as a completely integrable system once we 
obtain four independent integrals of motion through this procedure. 



3. Connection between integrating factors and nature of 

equations 

We note that equations (2.7)-(2.13) constitute an overdetermined system for the 

four unknowns, namely S, U, R and K. Among these equations, the first four equa- 
tions, (2.7)-(2.10), constitute the evolution equations for the variables S, U,R and 
K. We mention here that combining equations (2.7)-(2.10) one can get the following 
two identities, namely 

D[RS] = -{R<l>i^ + K(j)2^), D[KU] = -{R<i>iy + K<l>2y). (3.1) 

Any solution {S,U,R,K) which satisfies the equations (2.7)-(2.10) also satisfies 

the equation (3.1). Once the functions S, U, R and K are fixed then the rest 
of the problem is to verify whether these functions satisfy the 'extra determining 
equations', that is (2.11)-(2.13), or not. If these functions satisfy the extra deter- 
mining equations then they form a compatible set of solution and one can proceed 
to construct the associated integral of motion from (2.14). On the other hand if 
the functions do not satisfy the extra determining equations then one has to look 
for alternate ways to obtain compatible solutions. In fact, in practice, one often 
meets the case in which certain solution(s) which satisfies(satisfy) the evolution- 
ary determining equations (2.7)-(2.10) does(do) not satisfy the extra determining 
equations. More specifically, for a class of problems one often gets one or two or 
even three sets of {S, U, R, K) by solving (2.7)-(2.10) which does(do) satisfy the rest 
and another (other) set(s) does(do) not satisfy the later equations. In this situation 
we find an interesting fact that one can use the intcgral(s) derived from the set(s) 
which satisfies (satisfy) all the equations (2.7)-(2.13) and deduce the other compat- 
ible solution(s) {S,U,R,K) (definition of R and K follows). For example, let the 
set (5*4, U4, i?4, K4) be a solution of the evolution equations (2.7)-(2.10) and not of 
the extra determining equations (2.11)-(2.13). After analyzing several examples we 
find that one can make the set (^4, U4, R4, K4) compatible by modifying the form 
of R4 and K4 as 

R = F{t,x,y,x,y)R, K = G{t,x,y,x,y)K, (3.2) 

where F and G are functions to be determined. However, after making these modi- 
fications, ^4 and U4 remain the same. A motivation to do this type of modification 
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came from our earlier work on the applicability of the PS procedure to scalar sec- 
ond order ODEs (Chandrasekar et al. 2005a) . However, unlike the scalar case, now 
wc have to choose two functions, namely F and G, appropriately, such that the 
compatible forms of R and K can be fixed. Since we have to choose two functions 
one might ask the question whether these two functions are necessarily different or 
the same. In the following wc answer this question. 

Since R and K should satisfy equations (2.9) and (2.10), wc have 

D[F]R + FD[R] = -{FR(j)ii + GK(j)2± + FRS), (3.3) 
D[G]K + GD[K] = -{FR(t)iy + GK(j)2y + GKU). (3.4) 

Substituting the expressions for D[R] and D[K] (vide (2.9) and (2.10)) in (3.3) and 
(3.4) and simplifying the resultant equations we find that 

D[F]R={F-G)Kc^2±, D[G]K = {G - F)R^,y. (3.5) 

On the other hand substituting the modified forms (3.2) in the integrability condi- 
tions (3.1), wc get 

D[F]RS = {F- G)K(j)2x, D[G]KU = (G - F)R(j)iy. (3.6) 

Combining equations (3.5) and (3.6) we obtain 

{F - G)K{Scj>2i - 02.) =0, (G - F)R{UcPiy - cp^y) = 0. (3.7) 

Prom equation (3.7) we can conclude that either F ^ G or F = G. In the first case, 
we have uncoupled equations, that is, 

F ^ G ^ ct>2± = (l>iy = <^2x = (t>iy = {R = FR, k = GK) (3.8) 

and in the second case we have coupled equations, that is, 

F = G^<t>2±,(^ivA2x.4>iv^^ {R = FR, k = FK). (3.9) 

Note that in the case oi F ^ G, the other possible solution S = (f)2x/<t'2x, U = 
(l>iy/4>ii), when used in (2.7) and (2.8) lead to inconsistencies and so this choice is 
not considered. Further, the case R = and K = leads to the trivial solution. The 
above analysis clearly shows that for the uncoupled equations one should choose 
F and G as different and for the coupled equations one should choose them as the 
same. 

Finally, we mention the important point that the functions F and G are nothing 
but functions of integrals of motion. To show this, substituting back the relations 
(3.8) and (3.9) in (3.5) and (3.6), respectively, we get 

D[F]R = 0, D[G]K = and D[F]RS = 0, D[G]KU = 0, ^ D[F] = = D[G]. 

(3.10) 

Thus, irrespective of whether the given equations arc coupled ones or uncoupled 
ones the functions F and G can always be taken as integrals of motion or functions 
of them and the problem now is how to determine them. In the following, we discuss 
this in detail. 
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4. Case 1: F ^ G (Uncoupled Equations) 

(a) Theory 

In this case, we have the fohowing form of equation of motion, namely 

X = (l>i{t,x,x), y = <j)2{t,y,y). (4.1) 

As a consequence the evolution equations for the functions S, U, R and K (vide 
equations (2.7)-(2.10)) are simplified to the forms 

D[S] = -d^i^ + S(j>ii + S\ (4.2) 

D[U] = -(t>2v + U<t)2y + U\ (4.3) 

D[R] = -RS-R<Pi±, (4.4) 
D[K] = -KU-K^2y- (4.5) 

Solving equations (4.2) and (4.3) one can get explicit forms of S and U. Substitut- 
ing the known forms of S and U into the equations (4.4) and (4.5) and integrating 
the resultant equations one can fix the forms of R and K. From the known forms of 
{S, U,R,K) one can fix the integrals of motion from (2.14). One can also consider 
equation (4.1) as two independent second order ODEs and solve the equations in- 
dependently by adopting the procedure given in our previous paper (Chandrasekar 
et al. 2005a). To avoid repetition we do not discuss the uncoupled case further here. 



5. Case 2: F — G (Coupled equations) 

(a) Theory 

Now we focus on the case F = G. To obtain the forms (5, U, R, K) one needs 
to solve the compatibility conditions (2.7) - (2.10). For this purpose, we rewrite 
equations (2.7)-(2.10) in terms of two variables alone, namely R and K, by elimi- 
nating S and U, and analyse the resulting coupled second order partial differential 
equations (PDEs) and obtain expressions for R and K. From the latter we deduce 
the forms of S and U through the relations (2.9) and (2.10). 

To rewrite equations (2.7)-(2.10) interms of R and K let us take a total derivative 
of equations (2.9) and (2.10). Doing so we get 

D^[R] = -D[R(t>ii + K(t>2i + RS], D^[K] = -D[K(j)2y + R<i)iy + KU]. (5.1) 

Using the identities (3.1), equation (5.1) can be rewritten in a coupled form for R 
and K as 

D'^[R] + D[R(l)^i + K(l)2i] = R(t>ia^ + K^2^, (5.2) 

D^[K] + D[K(l>2y + Rct>ly] = Rct>ly + K(f>2y. (5.3) 

One can note that the above determining equations (5.2) and (5.3) form a sys- 
tem of linear PDEs. To solve equations (5.2) and (5.3) one may assume a specific 
ansatz either polynomial or rational in x and y for R and K and substituting 
the known expressions of (pi and ^2 and their derivatives into (5.2) and (5.3) 
and solving them one can obtain expressions for the integrating factors R and 
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K. Once R and K are known then the functions (5, U) can be fixed through the 
relation (2.9)-(2.10). Knowing S,U,R and K, one has to make sure that this set 
(5, U, R, K) does also satisfy the remaining compatibility conditions (2.11) - (2.13). 
The set {S,U,R,K) which satisfies all the equations (2.7) - (2.13) is then the ac- 
ceptable solution and one can then determine the associated integral / using the 
relation (2.14). For complete integrability we require four independent compatible 
sets {Si,Ui,Ri,Ki), z= 1,2,3,4. 

As discussed in §3, suppose the sets {Si,Ui,Ri,Ki), i = 1,2,3, are found to 
satisfy the equations (2.7) - (2.13) and the fourth set {84,1/4, R4,K4) docs not 
satisfy equations (2.11) - (2.13). In this case to identify the correct form of R4 and 
K4 one may assume that R4 = F{Ii,I'2,Iz)Ri and K4 = F{Ii,l2,Ii)K4, where 
F{Ii, I2, 13) is a function of the integrals /i, I2 and I^. To determine the explicit 
form of F(/i, J2, /a) we proceed as follows. Substituting 

R4 = F{h,l2, 13)R4, K4 = F{h,h, h)K4 (5.4) 
into equations (2.11)-(2.13), we get the following set of six relations: 

(aif{ + 5iF^ + cif^) 

{asFj + hFJ^ + csF^) 
{a^F[ + + csf^) 

where 



{a2F[ + b2F!, + C2F!,) 
d2 

{ajFj + bjFf. + CjF^) 

d4 

{aeFj + beF^ + cqF^) 
de 



F, 



(5.5) 



ai = R{Iix - SIi±), bi ^ R{l2x - Shx), ci ^ R{l3x - Shi), 

di = {SR± + RS±~Rx), a2 = K{hy-Uhy), b2 = K {hy - U hy) , 

C2 = K{hy-Uhy), d2={UKy+KUy-Ky), ^ {RI ly ' U K h^) , 

b3 = {Rl2y-UKl2i), C3 = {Rhy-UKhi), d3 = {UKi + KUi-Ry), 

a4 = {Khx-SRhy), b4 = {Kl2x-SRl2y), 04 = {Kh. - SRhy), 
d4 = iSRy + RSy-Kx), a, = (Rhy-Khi), b5 = {Rl2y-Kl2±), 
C5 = {Rhy-Kh±), d5^{Ki-Ry), a^^{SRhy-UKhx), 
66 = {SRhy - UKhx), C6 = {SRhy - UKhx), 
de = {UKx + KUx - SRy - RSy) 

are all known functions of t, x, y, x and y and F- = dF/dli. 

Equation (5.5) represents an over determined system of equations for the un- 
known F. A simple way to solve this equation is to uncouple it for F^, (= ^) 
i = 1,2,3, and solve the resultant equations. For example, eliminating F2 and ^3 
from the first three relations in equation (5.5) we obtain an equation for F{ in the 
form 

^ (rfiC2 - Cid2){biC3 - 63C1) - (diC3 - Cid3){biC2 - 62C1) 

F (aiC2 - cia2)(6iC3 - 63C1) - (aiCs - cia3)(6iC2 - 62C1) " 
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On the other hand ehminating F{ and from equation (5.5) (again from the first 
three relations) we arrive at the following equations for F2 in the form 



F2 _ (diC2 - cirf2)(aiC3 - aaci) - (diCs - cid3)(aiC2 - a2Ci) 
F (61C2 - ci62)(aiC3 - asci) - (61C3 - ci63)(aiC2 - a2Ci) ' 

In the similar way one obtains the following equation for F3 

-F3 _ (rfib2 - &i(^2)(ai?>3 - aah) - {diba - b-id3){aib2 - a2?>i) 
F (C162 - biC2){aib3 - asbi) - {cibs - 6iC3)(ai62 - 0261) " 



(5.7) 



(5.8 



One can easily check that the combination of other relations in equation (5.5) 
along with the forms (5.6)-(5.8) gives rise to relations which are effectively nothing 
but the constraint equations (2.11)-(2.13) and so no new constraint actually arises. 
Consequently, equations (5.6)-(5.8) can be written as 

^=gi{h,h,h)F, ^=g2{h,l2,h)F and ^ = gs{h,l2,h)F, (5.9) 
oh 0I2 0I3 

where (?j's, i = 1,2,3, are functions of Ji, I2 and Iz- Now solving equations (5.9) 
one can obtain the explicit form of F{Ii,l2,Iz). Once F is known we can obtain 
the complete solution i?4 and K/^ from which, along with ^4 and Ui, the fourth 
integral li can be constructed using the expression (2.14). 

Finally, we note that in some cases one can meet the situation that the sets 
[Si, Ui, Ri, Ki), i — 1,2, alone arc found to satisfy the equations (2.7) - (2.13) and 
the third set (^a, Uz,Rz,Kz) (as well as the fourth set) does not satisfy equations 
(2.11) - (2.13). In this case F may be a function of the integrals I\ and I2 which can 
be derived from the sets {Si, Ui,Ri,Ki), i = 1,2. Wc need to find the explicit form 
of F(/i, 12) in order to obtain the compatible solution (^a, UsjRsjKs). To recover 
the complete form of R3 and K3 one may assume that R3 = F(7i,/2)i?a and 

= F{Ii, /2)-ft'3 and proceed as before and obtain the determining equations for 
F. Since is a function of Ii and I2 alone one essentially gets the same form (5.5) 
but without the factor F3 (= §j^)- Since Cj's i = 1, ...6, are also exclusive functions 
of /a and their derivatives they do not appear in the determining equations. Solving 
the resultant determining equations one can fix the form of F which in turn provides 
us .Ra and K3 from which one can constract the third integral 7a for the given 
problem. Now from the known forms /i,/2 and /a, one can proceed as before to 
obtain the fourth compatible set {84,1/4, R4, K4) and thereby obtain the fourth 
integral I4 also. 

Similarly, if the set {Si, Ui, Ri, Ki) alone is found to satisfy the equations (2.7) 
- (2.13) and the second set {82,1^2, R2, K2) does not satisfy equations (2.11) - 
(2.13) then, in this case the determining equations for F takes the form (5.5) with 
F2 = F^ = hi = Ci = 0, i = \, ...6. Following the procedure mentioned above one 
can first derive the second integral. From a knowledge of (/i,/2), one can again 
follow the above procedures and construct 7a and I 4. In the following we illustrate 
the theory with specific examples. 

(6) Example 2: Coupled modified Emden equation 

The modified Emden equation (MEE) and its variants arise in different branches 
of physics (Erwin et al. 1984; Dixon & Tuszynski 1990) and considerable attention 
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has been paid recently to explore the mathematical and geometrical properties of 
MEE and its variant equations (see for example, Mahomed & Leach 1985; Steeb 
1993; Chandrasekar et al. 20056; Euler et al. 2007). The equation reads 

x + 3kxx + kx^ + Xx = 0, (5.10) 

where k and A are arbitrary parameters, which can be completely integrated in 
terms of elementary functions (Chandrasekar et al. 2005&). It exhibits some un- 
usual features such as amplitude independent frequency, nonstandard Lagrangian 
and time independent Hamiltonian forms and transition from periodic to front-like 
solutions as the sign of the parameter changes (Chandrasekar et al. 20056) . 

We now consider a coupled two dimensional version of (5.10) and demonstrate 
how it can be integrated out by following the theory given in the previous section. 
The coupled system of equations reads 

X + 2{kix + k2y)x + (fcii + k2y)x + {kix + k2y)'^x + Xix — 0, 

y + 2{kix + k2y)y + {kix + k2y)y + (fcix + k2y)'^y + X2y = 0, (5.11) 

where ki and Xi, i = 1,2, are arbitrary parameters. 

To determine the integrating factors we seek a rational form of ansatz for R and 
K in the form (suggested by the equation (5.10)) 

a, + a2x + a,y ^ b,+b,i + b,y ^ 

(a4 + 05a; + aeyji (04 + a^x + agyy 

where q and r are arbitrary numbers (to be fixed) and a^'s, i = 1,2. ..6 and 6i's, 
i = 1,2, 3, are arbitrary functions of t, x and y. Substituting (5.12) into (5.2) and 
(5.3) and equating the coefficients of different powers of x and y to zero we get a 
set of linear PDEs for the variables a;'s, i = 1,2. ..6, and 6j's, j = 1,2,3. Solving 
the resultant equations we find the following sets of particular solutions for R and 
K: 



Ri 



'\Xi{k2y + fc|j/^ + A2) + kix{k2Xiy - AaV^)) 
g{x,y,x,yy 



_ e-'^^Xik2{x + kix'^ + x{y/^ + k2y)) 

9[x,y,x,y)^ 



"^'* A2/C1 j y + k2y^ + y{V^ + kix)) 
g{x,y,x,y)^ 



i?2 = 



^ ^ e-v^«(A2(fcii- + k'fx^ + Al) + k2y{X2kix - XiV^)) 
^ 9{x,y,x,yY 

Rz = -, i^3 = 0, (5.15) 

g{x,y,x,yY 

i?4 = 0, = -, (5.16) 

g{x,y,x,yY 

where g{x, y, x, y) = {kiX2X + fc2Ai?) + k1X2X^ + fc|Aiy^ + fcifc2(Ai + X2)xy + A1A2). 
Now substituting the above forms of i?i's and /sTj's, i = 1,2,3,4, into (2.9) and 
(2.10), we can obtain the corresponding Si's and Ui's, i = 1,2,3,4. As a conse- 
quence now we have four sets of independent solutions {Si, Ui, Ri, Ki), i = 1,2, 3, 4, 
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for the equations (2.7)-(2.10). Now we check the compatibihty of these solutions 
with the remaining equations (2.11)-(2.13). We find that only the first two sets 

{Si,Ui, R.i, Ki), i = 1,2, satisfy the extra constraints (2.11)-(2.13), and become 
compatible solutions. Substituting their forms separately into equation (2.14) and 
evaluating the integrals we get 



J. _ e-v^«(i; + (kix + k2y)x + V^X[x) 
^ 9{x,y,x,y) 
J. ^ e-^^\y + {kix + k'2y)y + V^y) 
^ 9ix,y,x,y) 



(5.17) 



However, the sets {Si,Ui, Ri, KiYs, « = 3,4, do not satisfy the extra constraints 
(2.11)-(2.13) which means that the form of R3 in the third set and K4 in the fourth 
set may not be in the 'complete form' (since K3 = R4 = 0) but might only be a 
factor of the complete form. 

To deduce the compatible set (6*3, U3, R3, K3), let us substitute R3 = F{Ii,l2) 
R3 into equations (5.6) - (5.8). As a result we get 

1 dF 

-hFi + F = 0, F^=0, f^=0, {Fi = —,i = l,2,3). (5.18) 

Upon integrating (5.18), we get F = l/lf, (the integration constants are set to zero 
for simplicity) which fixes the form of R3 as 

^ {x + {kix + k2y)x + \/-\ixy ^ ^ 

Now one can easily check that the set [S^.Uz, R3, K3) is a compatible solution for 
the full set of determining equations (2.7)-(2.13) which in turn provides I3 through 
the relation (2.14) as 

^ e^'^Hi + (hx + k2y)x - V^x) 
(.i + {kix + k2y)x + ^-Xix) 

Finally, in the fourth set we take K4 = F{Ii, I2, l3)K4 and substitute it into 
equations (5.6)-(5.8), to obtain 

Fi = 0, ^I2F^ + F = 0, Fi = 0. (5.21) 

Upon integrating (5.21), we get F = l/if which fixes the form of K4 as 

= - ^ (5.22) 

{y + {kix + k2y)y + \/-X2yY ' 

Now the set {S4, U4, R4, K4) is a compatible solution for the equations (2.7)-(2.13) 
which in turn provides I4 through the relation (2.14) as 

, _ e2v^*(2/ + (fcia; + k2y)y - \f^y) ^3) 



{y + {kxx + k2y)y + \f^y) 
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Using the explicit forms of the integrals, Ii, I2, h and I4, the general solution to 
equation (5.11) can be deduced directly as 



x{t) - 



where Ai = \/—Xi + V— A2, A2 = \/— Ai — \/— A2, A = /i/sfci A2 and 72 = I2lik2\i- 
Finally wc note that for the case Aj > 0, i = 1,2, the above general solution 
becomes a periodic one, that is 

l-4^cosM + 5i)-^cosM + .52)' ^ 

Bsm.{ijj2t + 52) , . 

~ 1 - 4^ cos(a;ii + <5i) - cos(a;2t + ^2) ' ^ ^ ^ 

where uj = y/X], j = 1,2, A = ujioj^he-'^^ B = W2Wl/2e"'^^ h = e'"^'^' and 
I4 = e~^'*^ are arbitrary constants and {Aki/ui + 5^2/^2) < 1- More details on 
the classical dynamics and characteristic features of this system will be published 
elsewhere. 

(c) Example 3: Coupled Mathews - Lakshmanan oscillator 

Let us consider the coupled Mathews - Lakshmanan oscillator discussed by 
Carinena et al. (2004a;2004fe;2007), namely 

\{x^ -\- + X(xy — yx)'^)x — a^x ^ 

^ = {VT>^) = '^^' 



where r = + and A and a are arbitrary parameters, Equation (5.27) is 
a two dimensional generalization of the one dimensional non-polynomial oscillator 
introduced by Mathews & Lakshmanan (1974), 

(1 + \x^)x - {Xx^ - o?)x = 0. (5.28) 

Now wc use the following ansatz for K and K to explore the integrating factors 
for the equation (5.27) in the form (suggested by the one dimensional oscillator 

(5.28) ) 

2 2 3 2 2 3 

i? = ai + a2X + osy + a^x + a^xy + oey + a-^x + a^x y + agy x + awy , 
K = 61 + 62a; + bsy + 642;^ + 652;^ + h&y^ + 670;^ + h^x^y + b^y^x + b^y^, 

(5.29) 

where o,'s and biS i = 1,2. ..10, are arbitrary functions of t, x and y. Substituting 

(5.29) into (5.2) and (5.3) and equating the coefficients of different powers of x 
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and y to zero and solving the resultant PDEs as before one obtains the following 
nontrivial solutions for R and K: 



Ri = y, Ki = -X, 

2A((1 + Xy^)x - Xxyy) 



Rl 

R3 = x + y + 
K3 = x + y- 
Ra ~^ 



2A((1 + \x'^)y - Xxyx) 



1 + Ar2 ' 1 + Ar2 

X{x + yfiXxyy - (1 + W)x) 

rTv2 ' 

A(a; + yY{{l + Aa;^)^ - Xxyx) 



1 + Ar2 



1 + Ar2 



tX[ {x + yf 



{x + yY(bi 



{1 + Xy^)x - Xxyy 



-{x + y){a^ + A(l + Xy{x + y))xy - A(l + Xx{x + y))f) 



-1 



1 + Ar2 



:.\2 



tXl (x + y) 



{x + yf(j)2 



(1 + Xx^)y — Xxyx 



-{x + y){a^ - Xy{l + X{x + y))x^ + A(l + Xx{x + y))xy) 



(5.30) 
(5.31) 

(5.32) 
(5.33) 
(5.34) 



Now substituting the above forms of RiS and KiS, i = 1, 2, 3, 4, into the equa- 
tion (2.9) and (2.10), we get the corresponding forms of SiS and t/^'s, i = 1, 2, 3, 4. 
As a consequence now we have four sets of independent solutions for the equation 
(2.7)-(2.10). Out of these we find that only the first three sets {i = 1,2,3) satisfy 
the extra constraints (2.11)-(2.13) and become compatible solutions. Substituting 
these forms separately into equation (2.14) and evaluating the integrals we arrive 
at 



Ji = {yx - xy), 



(g^ - A(a:^ +y'^ + X{yx - xyf')) 
1 + Ar2 ' 



J =(i; I yf I jx + yfia^-Mi^ + y'^ + Hyx-xy)'^)) 

1 + Ar^ 



(5.35) 
(5.36) 



However, the set {84,114, R4, K4) does not satisfy the extra constraints (2.11)- 
(2.13) which means that the forms of R4 and K4 are incomplete. So we assume that 
R4 and K4 are in the forms (5.4). As a consequence equations (5.6) - (5.8) lead us 
to the following determining equations for the unknown F, namely 

dF 

Fi = 0, 2l2F^ + F = 0, hF^ + F = 0, {Fi = —, i = 1,2,3). (5.37) 

Upon integrating (5.37) we get, F = , (the integration constants are set to 
zero for simplicity) which fixes the form of R4 and K4 as 



i?4 



R4 



(5.38) 



where I1 and I3 are given in equations (5.35) and (5.36). Now one can easily check 
that the set (^4, Ua, Ra, Ka) is a compatible solution for the equations (2.7)-(2.13) 
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which in turn provides I4 through the relation (2.14) in the form 
h = -1 



i2(a2 - A(a;2 + y"^ + X{yx - xy^)) 



1 + Ar2 



+tan ^ 



(x + 2/)2(l + Ar2) 



. (5.39) 



Note that the above integral is given for the first time in the literature. 

Using the explicit forms of the integrals /i, /2, and 14, the solution to 
equation (5.27) can be deduced directly as 



x{t)-. 



la^ - h 


- 2Xh - 


A2J2 


2A/2 


/a2 - h 


- 2A/3 - 


A2/2 



2A72 



sin(v/i2t + <5), (5.40) 



wnereo-i4 2^^^* L /3(72_„2+2a/3+a=/?) J 



(c?) Example: 4 Known Two-dimensional Integrable Systems 

We have also studied the integrability properties of certain well known two- 
dimensional nonlinear Hamiltonian systems, namely Henon-Heiles system (Ramani 
et al. 1989; Lakshmanan & Sahadevan 1993; Lakshmanan & Rajasekar 2003) and 
generalized van der Waals potential equation (Ganesan & Lakshmanan 1990; Lak- 
shmanan & Rajasekar 2003) through the modified PS method. Our analysis shows 
that these systems do admit only the known integrable cases. In the following we 
present the salient features of our analysis. 

(i) Henon-Heiles system 

Let us consider the generalized Henon-Heiles system (Henon & Heiles 1964) 

X = -{Ax + 2axy) y = -{By + ax"^ - jSy'^) (5.41) 

where A,B,a and [3 arc arbitrary parameters. We obtain the integrating factors 
Ri = X and Ki = y for arbitrary values of the parameters using the ansatz (5.29). 
Substituting the forms R\ and Ki into the equation (2.9) and (2.10), we get the 
corresponding forms of Si and Ui. Further, we find that this set also satisfies the 
extra constraints (2.11)-(2.13) and become a compatible solution. Substituting these 
forms into the equation (2.14) and evaluating the integrals we get 

h = \{x' + f+ Ax" + By'') + ax'y - ^y\ (5.42) 

which is nothing but the total energy. Further, we find that only for the following 
specific parametric choices {1) A = B, a = —(3, (ii) A, B arbitrary, (3 = —6a and 
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(iii) B = IQA, (3 = —16a, one can get additional set of integrating factors, that is, 
(i) R2 = -y, K2 = -X, 

{ii) R2 = 2xy -xy- (^—^^^ 2;, K2 = -xx, 

{iii) R2 = a{xy — 6yx)x^ — 3{x^ + Ax'^)x, K2 = axx^, 

respectively. 

Finding the corresponding forms of S and U from (2.9) and (2.10) and substi- 
tuting the above forms of R and K along with the corresponding S and U into 
(2.14) and evaluating the integrals we obtain the following second integrals for the 
above three parametric choices 

(i) I2 = xy + {Ay + a{y'^ + ^x'^))x, 

{ii) h = 4:{xy - xy)x + {AAy + ax^ + Aay'^)x'^ + — —{x^ + Ax^), 
{iii) I2 = 9(a;^ + Aa;^)^ + 12axx'^{3yx - xy) - 2a'^x^{6y'^ + x'^) - 12aAx^y, 
respectively. 

(ii) Generalized van der Waals potential 

The generalized van der Waals potential equation in two dimensions is given by 

^ = -(270;+^), y = -{2j(3^y+^), r=(a;2+y2)l (5.43) 
where 7 and /3 are arbitrary parameters which can be derived from the Hamiltonian 

H = l{x^ + f)+j{x'+PV)--- (5.44) 

Repeating the above procedure we find that the system (5.43) admits integrating 
factors 



(5.45) 



(i) 


R2 


= -y, 


K2 = X, 




= 1 


{ii) 


R2 


= xy- 2yx, 


K2 = XX, 




= 4, 


{iii) 


R2 


= yy, 


K2 = yx- 2xy, 


0' 


1 

~ A' 



and the corresponding second integrals of motion are (with 7 arbitrary) 



(0 h 


= {yx - 


xy), 


{tl) h 


= {xy- 


yx)x 


{iii) h 


= {yx- 


xy)y 



r 

,,2 



£-2|£, = (5.46) 

The integrals Ii and I2 ensure the integrability of the above systems for the 
respective parametric choices through Liouville theorem. 
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6. Method of deriving general solution 

In the previous sections we derived the integrals of motion by deducing the functions 
S, U, R and K. However, in some situations it is very difficult to solve the determin- 
ing equations and obtain explicit forms for these functions. Now wc assume a situa- 
tion in which we are able to derive only two sets of solutions {Si,Ui,Ri, Ki), i = 1,2, 
and their associated integrals of motion and the remaining two sets of functions are 
difficult to derive from the determining equations. Under this situation, one may be 
able to deduce, in many cases, the remaining two integrals from the known integrals 
themselves without worrying further about the functions {Si,Ui, Ri, Ki), i = 3,4. 
The underlying idea is the following. 

From the known integrals we deduce a set of co-ordinate transformations and 
using these transformations we rewrite the original ODEs as a set of two first order 
ODEs. By integrating the latter wc obtain the remaining two integration constants. 
The solution to the original problem can then be obtained just by inverting the 
variables. In the following, we first describe the theory and then illustrate the ideas 
with an example. 

We start the procedure with two known integrals of motion, namely 

h=T{t,x,y,x,y), and h = Q{t,x,y,x,y). (6.1) 

Now we split the functional form of the integrals Ii and I2 into two terms such 
that one involves all the variables {t,x,y,x,y) while the other excludes x and y, 
that is, 

h = Fi{t,x,y,x,y) + F2{t,x,y), h = F^{t,x,y,x,y) + Fi{t,x,y). (6.2) 
Let us split the functions Fx and ^3 further in terms of two functions, that is, 

^-^^ ( c.(t,.!,,.,^) l^-(^'-'^)) +^^(^^(*'-'^))' 

Now identifying a set of new variables in terms of functions G\, G2, G3 and G4 

as 

wi = Gi{t,x,y), zi= G2{t' ,x,y,x,y)dt' , 

J o 

W2 = G3{t,x,y), Z2= Gi{t' ,x,y,x,y)dt' , (6.4) 

J o 

equation (6.3) can be rewritten in the form 

/l=i^l(^)+i^2(«;i), I2=F3(^^+F,{W2). (6.5) 

In other words, 

F,{^^=h-F2{w,), F,[^^=h-F2{w2). (6.6) 
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Rewriting equation (6.6) one arrives at the decoupled equations 

^ =/,(.,), 'i^=hM, (6.7) 

azi az2 

which can in principle be integrated and two additional integration constants iden- 
tified. Now rewriting the solution in terms of the original variables one obtains the 
general solution for the given coupled second order ODE. Thus the new variables 
Wi, W2, zi and Z2 correspond to transformations which effectively decouple the 
original coupled nonlinear second order ODEs (2.1), provided such transformation 
variables can be identified. In the following we illustrate the procedure with an 
example. 

(a) Example: 5 Force-free coupled Duffing-van der Pol oscillators 
Let us consider the force-free coupled DufHng-van der Pol oscillators of the form 

x + {a + (3{x + yY)x - ujix + 5ix^ + 71 xy^ + Aix^y = 0, 
y+{a + l3{x + yf)y - uj^y + 52y^ + J2x'^y + X2xy'^ = 0, (6.8) 

where a, (3, uji, Si, -fi and ~ 1,2, are arbitrary parameters. Though the general 
equation (6.8) has not been discussed in detail in the literature, special cases of the 
above coupled Duffing-van der Pol equation have been used to represent physical 
and biological systems (Linkens 1974; Datardina & Linkens 1978; Kawahara 1980; 
Rajasekar & Murali 2004). 

We seek a rational form of ansatz for R and K in the form given by (5.12). 
Upon solving the equations (5.2)-(5.3) with this ansatz we find that three non- 
trivial solutions exist for the choice 

LOi = — (5i =7i = — and Ai = — , I = 1,2. (6.9) 

The corresponding forms of R and K turn out to be 

[U,,K,) [e ,e ), [H2,K2) ( (y+ay)' (y+|y)2^ 

Now substituting the form of Ri's and Ki's, i = 1,2, 3, into the equations (2.9) and 
(2.10), we get the corresponding S'^'s and Ui's, i = 1,2, 3. As a consequence now we 
have three sets of independent solutions for the equations (2.7)-(2.10) which are all 
found to be compatible with (2.11)-(2.13). Substituting the forms {Si, Ui, Ri, i^j)'s, 
i = 1,2, 3, separately into equation (2.14) and evaluating the integrals we get 



Ii = {x + y+-{x + y) + -{x + yy)e^ 



^ {x + fx) ^ / (xij - yx) 



e2"*. (6.11) 
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To illustrate the theory given in the first part of this section we consider the 
first two integrals Ii and I2 given by equation (6.11), and rewrite them in the form 
(6.2) as 



h = 
h = 



{x + \x)e 



\oct 



(2/+fy)e4«* 

Now splitting the first term in I\ and /2 further in the form (6.3) as 



(6.12) 



d 



e^"*(^^[(a: + y)e^"*] j +|(x + y)^e 



iat 



(6.13) 



we identify the new dependent and independent variables from (6.13) by using the 
relations (6.4): 



Wi = [X + y)e* 



2 -iat 

Zl = e 

a 



W2 = xe'^ , 



Z2 = ye" 



(6.14) 



One can easily check that in the new variables equations (6.11) become uncoupled 
as 



pwi— — = 0, 



= 0. 



(6.15) 



d^zi dzi ' d'^Z2 

The integrals of (6.15) can be easily obtained. These are nothing but the inte- 
grals Ii and I2 given in (6.12) but in terms of the new variables (6.14) they read 
as 



dwi , /3 3 J dw2 
dzi 6 az2 
Solving equations (6.16), we obtain (Gradshteyn & Ryzhik 1980) 



(6.16) 



W2 = I2Z2 + h, 



1 



log 



(wi + of 
w\ — awi + 



+ VS tan-i 



Wi 



73 



2a — wi 



(6.17) 



where a = y^^i h and zq are the third and fourth integration constants. Rewrit- 
ing Wi and Zi, i = 1,2, in terms of the old variables, one can get the explicit solution 
for equation (6.8) for the parametric choice (6.9) in the form 

2 _^ a \1 f iix + y)e>' + a)^ 

--e -zo = x^- klog ' 
a ill I z 



-|-\/3 tan"^ 



(x + j/)2e2"* — a{x + y)e4"* + 
(a; + 2/)e3"*y3 ^ 



2a — {x + y)e^^ 



(6.18) 



x = l2y + he-^"\ 

from which one can obtain implicit relations between x and t and y and t, corre- 
sponding to the general solution. 
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7. Conclusion 

In this paper, as a first part of our investigations on tlie complete integrability and 
linearization of two coupled second order ODEs, we have focused our attention only 
on the intograbihty aspects. In particular, we have introduced a general method of 
finding integrable parameters, integrating factors, integrals of motion and general 
solution associated with a set of two coupled second order ODEs through the ex- 
tended PS procedure. The procedure can also be extended straightforwardly to 
analyse any number of coupled second order ODEs. The proposed method is sim- 
ple, straightforward and very useful to solve a class of coupled second order ODEs. 
We have illustrated the theory with potentially important examples. We have also 
introduced a novel idea of transforming coupled second order nonlinear ODEs into 
uncoupled second order ODEs. We have deduced the transformations from the first 
two integrals themselves. 

It is also of interest to study the problem of linearization of coupled nonlinear 
ODEs by transforming them into linear ODEs as in the case of single second or- 
der nonlinear ODEs. However, it is a more difficult and challenging problem than 
quadratures. The primary reasons are (1) First of all it is not known, in general, 
whether the given equation is linearizable or not since in the case of two degrees 
of freedom systems there can exist many types of linearizing transformations and 
it is not obvious which one will be successful. (2) What are the possible transfor- 
mations that could exist in the case of two coupled second order ODEs is also not 
known in the literature. (3) Finally, there is no simple procedure that gives us the 
required transformation in a straightforward way. We would like to address all these 
questions in the follow up paper (Paper V). 

The work of MS forms part of a research project sponsored by National Board 
for Higher Mathematics, Government of India. The work of ML forms part of a 
Department of Science and Technology, Government of India sponsored research 
project and is supported by a DST Ramanna Fellowship. 
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